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ABSTRACT 

We present the results of a survey of the Coma Berenices open star cluster (Melotte 111), 
undertaken using proper motions from the USNO-B1.0 and photometry from the 2MASS 
Point Source catalogues. We have identified 60 new candidate members with masses in the 
range l.OO7M <M<O.269M0. For each we have estimated a membership probability by 
extracting control clusters from the proper motion vector diagram. All 60 are found to have 
greater than 60 per cent probability of being clusters more than doubling the number of known 
cluster members. The new luminosity function for the cluster peaks at bright magnitudes, 
but is rising at K« 12, indicating that it is likely lower mass members may exist. The mass 
function also supports this hypothesis. 
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1 INTRODUCTION 

During the last two decades there have been numerous deep surveys 
of young nearby open clusters focusing on the detection of very low 
mass stellar and substellar members (e.g. Jameson & Skillen 1989, 
Lodieu et al. 2005). Since these objects fade during their evolu- 
tion, in these environments they are comparatively luminous, plac- 
ing them comfortably within the reach of 2/4m class telescopes. 
Furthermore, numerical simulations suggest that in clusters with 
ages less than ~200Myrs, dynamical evolution should not yet have 
led to the evaporation of a large proportion of these members (e.g. 
de la Fuente Marcos & de la Fuente Marcos 2000). 

For these same reasons, there have been relatively few searchs 
of this type undertaken in older open clusters. Indeed, the few, deep, 
large area surveys of the Hyades, the most extensively studied clus- 
ter of its type and the closest to the Sun (46.3pc; Perryman et al. 
1998), have led to the identification of only a very small number 
of low mass stellar and substellar members (e.g. Reid 1993, Gizis 
et al. 1999, Moraux, priv. comm.). This finding is loosely in agree- 
ment with the predictions of N-body simulations which indicate 
that less than a fifth of the original population of substellar mem- 
bers remains tidally bound to a cluster at the age of the Hyades, 
625±50Myrs (Perryman et al. 1998). However, despite this rough 
consistency, until the very low mass populations of more open clus- 
ters of similar age to the Hyades have been studied it seems prema- 
ture to completely exclude other interpretations for the deficit of 
low mass members here e.g. differences in the initial mass func- 
tion. Furthermore, additional investigations of this nature may also 
be of use in refining N-body simulations. 

Therefore we have recently embarked on a survey of the Coma 
Berenices open star cluster (Melotte 111, RA = 12 23 00, Dec = 
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+26 00 00, J2000.0) to extend our knowledge of it's luminosity 
function towards the hydrogen burning limit. At first glance this 
seems a prudent choice of target. It is the second closest open clus- 
ter to the Sun. Hipparcos measurements place it at d=89.9±2.1pc 
(van Leeuwen 1999), in agreement with older ground based esti- 
mates (e.g. d=85.4±4.9pc, Nicolet 1981). Furthermore, foreground 
extinction along this line of sight is low, E(B-V)~0.006±0.013 
(Nicolet 1981). The metalicity of the cluster is relatively well con- 
strained. Spectroscopic examination of cluster members reveals it 
to be slightly metal poor with respect to the Sun. For example, 
Cayrel de Strobel (1990) determine [Fe/H] =-0.065±0.021 using a 
sample of eight F, G and K type associates, whereas Friel & Boes- 
gaard (1992) determine [Fe/H]=-0.052±0.047 from fourteen F and 
G dwarf members. While estimates of the age of Melotte 111 vary 
considerably from 300Myrs to lGyr (e.g. Tsvetkov, 1989), more 
recent determinations, based on fitting model isochrones to the ob- 
served cluster sequence, are bunched around 400-500Myrs (e.g. 
Bounatiro & Arimoto 1993, Odenkirchen 1998). Thus the Coma 
Berenices open cluster is probably marginally younger the Hyades. 

However, Melotte 111 is projected over a large area of sky 
(~100 sq. deg.) and contains considerably fewer bright stellar 
members than the Hyades. For example, Odenkirchen (1998) deter- 
mines the cluster tidal radius to be ~5-6pc but finds only 34 kine- 
matic members down to V=10.5 within a circular area of radius 5° 
centred on the cluster. He estimates the total mass of Melotte 111 
to lie in the range 30-90Mq , which can be compared to estimates 
of 3OO-46OM for the mass of the Hyades (e.g. Oort 1979, Reid 
1992). Additionally, the small proper motion (jU a =-11.21 ±0.26 
mas yr~\ /i<s=-9.16 ±0.15 mas yr _ ; van Leeuwen 1999) means 
that proper motion alone is not a suitable means by which to dis- 
criminate the members of Melotte 111 from the general field pop- 
ulation. Fortunately, the convergent point for the cluster is suffi- 
ciently distant at a = 6 40 31.2, S = -41 33 00 (J2000)(Masden et 
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al 2002), that we can expect all the cluster members to have essen- 
tially the same proper motion. 

In the first detailed survey of Melotte 111, Trumpler (1938) 
used proper motion, spectrophotometric and radial velocity mea- 
surements to identify 37 probable members, mp<10.5, in a cir- 
cular region of 7° diameter centered on the cluster. A significant 
additional number of fainter candidate members were identified 
by Artyukhina (1955) from a deep (mp<15) proper motion sur- 
vey of ~7 sq. degrees of the cluster. Argue & Kenworthy (1969) 
performed a photographic UBVI survey of a circular field, 3.3° in 
diameter, to a limiting depth of mp=15.5. They rejected all but 2 of 
her candidates with mp>ll but identified a further 2 faint objects 
with photometry and proper motion which they deemed to be con- 
sistent with cluster membership. Subsequently, De Luca & Weis 
(1981) obtained photoelectric photometry for 88 objects (V>11), 
drawn from these two studies. They concluded that only 4 stars, 
3 of which were listed in Argue & Kenworthy as probable mem- 
bers, had photometry and astrometry consistent with association to 
Melotte 111. 

More recently, Bounatiro (1993) has searched an area of 
6°x6° centered on Melotte 111, using the AGK3 catalogue, and 
has identified 17 new candidate members (mp<12). However, de- 
spite Randich et al. (1996) identifying 12 new potential low mass 
members (Vwll.5-16) from a ROSAT PSPC survey of the cluster, 
a detailed follow-up study has concluded that none of these are as- 
sociates of Melotte 111 (Garcia-Lopez et al. 2000). Odenkirchen 
et al. (1998) have used the Hipparcos and ACT catalogues to per- 
form a comprehensive kinematic and photometric study of 1200 sq. 
degrees around the cluster center complete to a depth of V^IO.5. 
They find a total of ~ 50 kinematic associates which appear to 
be distributed in a core-halo system. The core, which is elliptical 
in shape with a semi-major axis of 1.6°, is dominated by the more 
massive members while the halo contains proportionately more low 
mass associates. Odenkirchen et al. also find evidence of an extra- 
tidal 'moving group' located in front of the cluster in the context of 
its motion around the Galaxy. However, from a subsequent spectro- 
scopic study, Ford et al. (2001), concluded that approximately half 
of the moving group were not associated with the cluster. 

The ready availability of high quality astrometric and photo- 
metric survey catalogues (e.g. USNO-B1.0, 2MASS) have made a 
new deeper, wide area survey of the Coma Berenices open star clus- 
ter a tractable undertaking. In this paper we report on our efforts to 
use the USNO-B 1.0 and 2MASS Point Source Catalogues to search 
for further candidate low mass members of Melotte 111. We iden- 
tify 60 new candidates with proper motions and photometry consis- 
tent with cluster membership. State-of-the-art evolutionary models 
indicate some of these objects have masses of only M^O.269 M . 



2 THE PRESENT SURVEY 

In ideal cases, membership of nearby open clusters can be ascer- 
tained by either photometry or proper motion. In the former ap- 
proach, associates of a cluster are generally identified by insisting 
that they have magnitudes and colours which, within uncertainties, 
sit them on or close to an appropriate model isochrone (e.g. Bouvier 
et al. 1998). In the latter approach, as measurement errors usually 
dominate over the cluster velocity dispersion, selection of mem- 
bers simply requires that objects have astrometric motions consis- 
tent with known associates of a cluster. When combined, these two 
methods can yield candidates with a very high probability of cluster 
membership (e.g. Moraux et al. 2001) 
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Figure 1. Plot of the UCAC (Zacharias et al. 2004) proper motions for 4 
square degrees of the cluster centre. The thin line borders the cluster proper 
motion selection, and the thick line, the contols. The shallower UCAC cata- 
logue is used here for illustration as the USNO B 1 .0 catalogue only provides 
quantized values of 2 mas yr -1 . 



However, as discussed in §1, the proper motion of Melotte 
1 1 1 is comparatively small and not suitable alone for the discrima- 
tion of cluster members. For example, members of the Pleiades 
open cluster, which can be readily identified astrometrically, have 
proper motions of /io,=19.14±0.25 mas yr~ , ^=-45.25 ±0.19 
mas yr -1 (van Leeuwen 1999). Nevertheless, because of the large 
epoch difference of approximately 50 years between the two Palo- 
mar Sky Surveys, if we restrict the search to R<18.0, the proper 
motion measurements available in the USNO-B 1.0 catalogue are 
sufficiently accurate that they can be used in conjunction with suit- 
able photometry to identify candidates with a substantial probabil- 
ity of being members of Melotte 111. In principle the USNO-B 1.0 
catalogue provides B, R and I photometry but the significant uncer- 
tainties in the photographic magnitudes (~0.3 mags) severely limit 
its usefulness for this work (Monet et al. 2003). Therefore, to com- 
pliment the astrometry, we choose to use J, H and Ks photometry 
from the 2MASS Point Source Catalogue which has S/N<:10 down 
to J=15.8, H=15.1 and K s = 14.3 (Skrutskie et al. 1997). 

The survey was conducted as follows:- 

(i) A circular area of radius 4 degrees centred on the cluster was 
extracted from the USNO B1.0 catalogue. 

(ii) Stars were selected according to the criterion (p a - X) 2 
+ (p s - Y) 2 < 100, where X = -11.21 and Y = -9.16, i.e. to 
lie within 10 mas yr -1 of the Hipparcos determined value for 
the proper motion of Melotte 111. This procedure was initially 
repeated for X = -11.21, Y = +9.16 and X = +11.21, Y = -9.16, 
to obtain two control samples (see Figure 1). Two more control 
samples were later extracted from two further circular regions 
of USNO B1.0 data. These data had a similar galactic latitude 
to the cluster but were offset by 10° from the centre of Melotte 
111. Stars were selected from these latter regions by applying the 
first proper motion criterion above. The known members have 
errors, or velocity dispersions ammounting to about ± 2.0 km s _1 , 
(Odenkirchen et al 1998). This corresponds to a proper motion of 
± 4.8 mas yr _1 , and the USNO B1.0 catalogue astrometry errors, 
(see Tables 1 and 2) are small - typically less than ± 5 mas yr 
_1 for our stars. Note the USNO B1.0 proper motion errors are 
quantized in units of 1 mas yr~ 1 and a zero error thus indicates an 
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error of less than 0.5 mas yr . Thus if we selected a bounding 
circle of 10 mas yr, and the total quadratically added error is 7 mas 
yr~\ we have selected all stars to a completeness level of 1.4 a, 
which means our survey is complete to w 90%. 

(iii) Stars passing the above sets of criteria were cross- 
referenced against the 2MASS point source catalogue using a 
match radius of 2 arcseconds. 

(iv) Subsequently the sample of candidate cluster members 
and the four control samples were plotted in Ks, J-Ks colour- 
magnitude diagrams as illustrated in Figures 2a and 2b respectively. 

A cursory glance at the CMDs in Figure 2 reveals that our 
method appears to be rather successful in finding new associates of 
Melotte 111. An obvious cluster sequence can be seen extending 
to K«12 beyond which it is overwhelmed by field star contamina- 
tion. To perform a quantitative selection of candidate members we 
restrict ourselves to Ks<12. We use as a guide to the location of 
the cluster sequence the previously known members and a 400Myr 
NEXTGEN isochrone for solar metalicity (Baraffe et al. 1998), 
scaled to the Hipparcos distance determination for Melotte 111. 
In the magnitude range where they overlap, the previously known 
cluster members of the single star sequence, congregate around the 
theoretical isochrone (all previously known cluster members are 
listed in Table 1). Furthermore, the model isochrone appears to con- 
tinue to follow closely the excess of objects in Figure 2a, relative 
to Figure 2b, suggesting that it is relatively robust in this effective 
temperature regime. 

The location of the theoretical isochrone at J-Ks< 0.8 is in- 
sensitive to the uncertainties in the age of the cluster. The cluster 
would have to be much younger to significantly shift the isochrone. 
The bulk of the observed dispersion in the single star sequence here 
likely stems from the finite depth of the cluster (~0. 15 mags). Nev- 
ertheless, in this colour range we choose to select objects which lie 
no more than 0.3 magnitudes below the theoretical isochrone as this 
ensures that all previously known cluster members are recovered 
(filled squares in Figure 2a). As binary members can lie up to 0.75 
magnitudes above the single star sequence, only objects which also 
lie no more than 1.05 magnitudes above the theoretical sequence 
are deemed candidate members. 

Redward of J-Ks=0.8, the main sequence becomes very steep 
in this CMD. We have tried using various combinations of 2MASS 
and USNO-B1.0 photometry (e.g. colours such as R-K) to circum- 
vent this. However, as mentioned previously, the poorer quality of 
the photographic magnitudes provided by the USNO-B1.0 cata- 
logue results in a large amount of scatter in optical+IR CMDs ren- 
dering them of little use for this work. Nonetheless, the finite depth 
of the cluster and a small error in cluster distance determination 
have a negligible effect in this part of the Ks, J-Ks CMD. Based 
on previous experience gained from our investigations of the low 
mass members of the Pleiades and Praesepe open clusters, we esti- 
mate an uncertainty in the model J-K colour of ± 0.05 magnitudes. 
Hence for J-Ks>0.8, we have selected all objects which overlap 
a region 0.1 magnitudes wide in J-K, centered on the theoretical 
isochrone. 

To assess levels of contamination in our list of candidate mem- 
bers, we have imposed these same colour selection criteria on the 
control samples. The resulting sequences were divided into one 
magnitude bins (in Ks)for J-Ks>0.8, and in bins of 0.2 in J-Ks for 
J-Ks <0.8 and the number of objects in each for both the cluster and 
control samples counted. Subsequently, the membership probablil- 
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Figure 2. The CMD for the cluster (top) and a control sample with 
/i Q =+11.21 ^=-9.16 (bottom). Previously known members of the clus- 
ter are highlighted in the upper plot (solid squares). A 0.4Gyr NEXTGEN 
isochrone is overplotted from Baraffe et al (1998) (solid line). This was con- 
verted into the 2MASS system using the transforms of Carpenter(2001). 



ity for each candidate cluster member, ^membership, was estimated 
using equation 1, 



membership 



N i 
1 ' cms 



-^control 



iVcl 



(1) 



where Nd^t^ is the number of stars in a magnitude bin from the 
cluster sample and N contIol is the mean number of stars in the same 
magnitude range but drawn from the control samples. Our list of 
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Table 1. Name,coordinates,proper motion measurements, R,I,J,H and K magnitudes for the known cluster members as detailed in Bounatiro and Arimoto 
(1992) and Odenkirchen et al. (1998) for the area we surveyed. The masses were calculated by linearly interpolating the models of Girardi et al(2002) M^l 
M Q orBaraffeetal(1998),M<l M . 



Name RA Dec. fi a fig error/i Q error fi$ R I J H Ks Mass 

J2000.0 masyr- 1 M 



BD+26 2337 


12 22 30.31 


+25 50 46.1 


-12.0 


-10.0 


0.0 


0.0 


4.54 


4.28 


3.78 


3.40 


3.23 


2.709 


BD+28 2156 


12 51 41.92 


+27 32 26.5 


-10.0 


-10.0 


0.0 


0.0 


4.55 


4.20 


3.62 


3.36 


3.26 


2.709 


BD+28 2115 


12 26 24.06 


+27 16 05.6 


-14.0 


-10.0 


0.0 


0.0 


4.79 


4.65 


4.40 


4.23 


4.14 


2.685 


BD+24 2464 


12 29 27.04 


+24 06 32.1 


-18.0 


0.0 


0.0 


0.0 


5.25 


5.02 


4.86 


4.57 


4.54 


2.595 


BD+27 2134 


12 26 59.29 


+26 49 32.5 


-14.0 


-10.0 


0.0 


0.0 


4.91 


4.88 


4.79 


4.72 


4.64 


2.568 


BD+26 2344 


12 24 18.53 


+26 05 54.9 


-14.0 


-14.0 


0.0 


0.0 


5.11 


5.08 


4.93 


4.94 


4.89 


2.496 


BD+25 2517 


12 31 00.56 


+24 34 01.8 


-12.0 


-10.0 


0.0 


0.0 


5.41 


5.39 


5.29 


5.30 


5.26 


2.363 


BD+26 2354 


12 28 54.70 


+25 54 46.2 


-24.0 


-16.0 


0.0 


0.0 


5.24 


5.28 


5.22 


5.29 


5.28 


2.356 


HD 105805 


12 10 46.09 


+27 16 53.4 


-12.0 


-12.0 


0.0 


0.0 


5.92 


5.86 


5.65 


5.66 


5.60 


2.216 


BD+26 2345 


12 24 26.79 


+25 34 56.8 


-14.0 


-14.0 


0.0 


0.0 


6.59 


6.46 


5.84 


5.77 


5.73 


2.150 


BD+23 2448 


12 19 19.19 


+23 02 04.8 


-14.0 


-10.0 


0.0 


0.0 


6.14 


6.04 


5.94 


5.96 


5.90 


2.057 


BD+26 2326 


12 19 02.02 


+26 00 30.0 


-14.0 


-10.0 


0.0 


0.0 


6.36 


6.27 


6.08 


6.00 


5.98 


2.018 


BD+25 2523 


12 33 34.21 


+24 16 58.7 


-12.0 


-10.0 


0.0 


0.0 


6.19 


6.13 


6.03 


5.98 


5.98 


2.014 


BD+27 2138 


12 28 38.15 


+26 13 37.0 


-16.0 


-8.0 


0.0 


0.0 


6.40 


6.30 


6.13 


6.02 


5.99 


2.011 


BD+26 2343 


12 24 03.46 


+25 51 04.4 


-14.0 


-10.0 


0.0 


0.0 


6.59 


6.47 


6.17 


6.07 


6.05 


1.978 


BD+26 2353 


12 28 44.56 


+25 53 57.5 


-22.0 


-18.0 


0.0 


0.0 


6.52 


6.40 


6.16 


6.10 


6.05 


1.977 


BD+29 2280 


12 19 50.62 


+28 27 51.6 


-12.0 


-10.0 


0.0 


0.0 


6.52 


6.42 


6.20 


6.19 


6.13 


1.931 


BD+26 2352 


12 27 38.36 


+25 54 43.5 


-14.0 


-12.0 


0.0 


0.0 


6.57 


6.48 


6.28 


6.22 


6.22 


1.879 


BD+30 2287 


12 31 50.55 


+29 18 50.9 


-12.0 


-10.0 


0.0 


0.0 


7.37 


7.20 


6.84 


6.74 


6.65 


1.607 


BD+25 2495 


12 21 26.74 


+24 59 49.2 


-12.0 


-10.0 


0.0 


0.0 


7.23 


7.08 


6.79 


6.74 


6.66 


1.600 


BD+26 2347 


12 25 02.25 


+25 33 38.3 


-14.0 


-8.0 


0.0 


0.0 


7.83 


7.55 


7.05 


6.85 


6.76 


1.540 


BD+26 2323 


12 17 50.90 


+25 34 16.8 


-12.0 


-12.0 


0.0 


0.0 


7.66 


7.47 


7.08 


6.98 


6.92 


1.451 


BD+26 2321 


12 16 08.37 


+25 45 37.3 


-12.0 


-10.0 


0.0 


0.0 


7.87 


7.66 


7.23 


7.11 


7.03 


1.399 


BD+28 2087 


12 12 24.89 


+27 22 48.3 


-12.0 


-12.0 


0.0 


0.0 


7.85 


7.64 


7.27 


7.13 


7.08 


1.380 


BD+28 2095 


12 16 02.50 


+28 02 55.2 


-24.0 


-6.0 


0.0 


0.0 


8.03 


7.80 


7.41 


7.22 


7.20 


1.331 


BD+27 2129 


12 25 51.95 


+26 46 36.0 


-14.0 


-10.0 


0.0 


0.0 


8.09 


7.86 


7.41 


7.30 


7.20 


1.331 


BD+27 2122 


12 23 41.00 


+26 58 47.7 


-14.0 


-10.0 


0.0 


0.0 


8.10 


7.87 


7.46 


7.33 


7.25 


1.313 


BD+23 2447 


12 18 36.17 


+23 07 12.2 


-14.0 


-10.0 


0.0 


0.0 


8.39 


8.13 


7.63 


7.38 


7.30 


1.294 


BD+28 2109 


12 21 56.16 


+27 18 34.2 


-10.0 


-12.0 


0.0 


0.0 


8.22 


7.96 


7.56 


7.39 


7.32 


1.285 


HD 107685 


12 22 24.75 


+22 27 50.9 


-12.0 


-10.0 


0.0 


0.0 


8.23 


8.00 


7.60 


7.39 


7.38 


1.263 


BD+24 2457 


12 25 22.49 


+23 13 44.7 


-14.0 


-10.0 


0.0 


0.0 


8.30 


8.06 


7.64 


7.48 


7.39 


1.261 


BD+28 2125 


12 31 03.09 


+27 43 49.2 


-16.0 


-8.0 


0.0 


0.0 


8.27 


8.01 


7.61 


7.46 


7.40 


1.256 


BD+25 2488 


12 19 28.35 


+24 17 03.2 


-12.0 


-12.0 


0.0 


0.0 


8.69 


8.40 


7.86 


7.55 


7.49 


1.224 


HD 109483 


12 34 54.29 


+27 27 20.2 


-12.0 


-10.0 


0.0 


0.0 


8.67 


8.40 


7.89 


7.58 


7.51 


1.217 


BD+25 2486 


12 19 01.47 


+24 50 46.1 


-12.0 


-10.0 


0.0 


0.0 


8.53 


8.27 


7.83 


7.55 


7.53 


1.207 


BD+27 2130 


12 26 05.48 


+26 44 38.3 


-8.0 


-14.0 


0.0 


0.0 


9.44 


9.04 


8.13 


7.68 


7.58 


1.193 


HD 107399 


12 20 45.56 


+25 45 57.1 


-12.0 


-8.0 


0.0 


0.0 


8.68 


8.42 


7.97 


7.74 


7.65 


1.171 


BD+26 2340 


12 23 08.39 


+25 51 04.9 


-12.0 


-10.0 


0.0 


0.0 


8.80 


8.52 


8.02 


7.76 


7.68 


1.160 


BD+25 2511 


12 29 40.92 


+24 31 14.6 


-10.0 


-10.0 


0.0 


0.0 


9.26 


8.80 


8.20 


7.84 


7.72 


1.148 


BD+27 2121 


12 23 41.82 


+26 36 05.3 


-16.0 


-12.0 


0.0 


0.0 


8.85 


8.49 


8.13 


7.79 


7.73 


1.143 


BD+27 2117 


12 21 49.02 


+26 32 56.7 


-12.0 


-8.0 


0.0 


0.0 


8.89 


8.57 


8.21 


7.86 


7.85 


1.109 


TYC1991-1087-1 


12 27 48.29 


+28 11 39.8 


-12.0 


-10.0 


().() 


().() 


9.26 


8.94 


8.43 


8.05 


8.05 


1.051 


HD 105863 


12 11 07.38 


+25 59 24.6 


-12.0 


-10.0 


0.0 


0.0 


9.14 


8.84 


8.38 


8.11 


8.07 


1.043 


BD+30 2281 


12 29 30.02 


+29 30 45.8 


10.0 


0.0 


0.0 


0.0 


9.41 


9.12 


8.38 


8.16 


8.07 


1.043 


BD+36 2312 


12 28 21.11 


+28 02 25.9 


-14.0 


-12.0 


0.0 


0.0 


9.92 


9.54 


8.94 


8.47 


8.46 


0.915 



candidate associates of Melotte 1 1 1 is presented in Table 2, along 
with these estimates of membership probability. 



We note that there is a slight increase in the level of contami- 
nation in the range 0.55<J-Ks <0.7. A similar increase in the num- 
ber of field stars was seen by Adams et al.(2002) when studying 
the Praesepe open cluster, which, like Melotte 111, has a relatively 
high galactic latitude, b=38°. We believe this is due to K giants 
located in the thick disc. 



3 RESULTS 

Our survey of the Coma Berenices open cluster has recovered 45 
previously known members in total, 38 listed by Bounatiro and Ari- 
moto (1992) and 7 unearthed by Odenkirchen et al (1998). Further- 
more, it has identified 60 new candidate cluster members with mag- 
nitudes down to Ks=12. Beyond this magnitude, no statistically 
significant difference between the cluster and the control samples 
was found, as the contamination by field stars is too great. We be- 
lieve that our survey is reasonably complete to this limit; we tried 
proper motion search radii of 7 and 5 mas yr _1 , but in both cases 
found increasing numbers of likely members were excluded. Ex- 
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Table 2. Coordinates.proper motion measurements, R,I,J,H and K magnitudes. Mass calculated from linearly interpolating the NEXTGEN model and the 
absolute K magnitude. The probability of membership for each of our 60 possible members. 



RA Dec. 
J2000.0 


Met 


fis error/i Q 
mas yr _1 


error fig 


R 


I 


J 


H 


K S 


Mass 

Mq 


Membership 
probability 


12 24 17.15 


+24 19 28.4 


-18 


-14.0 


0.0 





9 44 


9 08 


8 42 


7 QS 

/ .y J 


7 90 


880 


64 


12 38 14.94 


+26 21 28.1 


-6 


-2.0 


0.0 


u.u 


Q S7 

y.J 1 


9 13 


8 56 


8 07 


8 00 


836 


64 


12 23 28.69 


+22 50 55.8 


-14 


-10.0 


0.0 





Q 77 

y. I l 


Q 91 


8 60 


8 09 


8 01 


815 


64 


12 31 04.78 


+24 15 45.4 


-8 


-4.0 


0.0 


u.u 


9 76 


Q 1 Q 

7.17 


8 81 


8 28 


8 20 


798 


79 


12 27 06.26 


+26 50 44.5 


-12 


-8.0 


0.0 





9 59 


9 31 


8 64 


8 33 


8 25 


1 007 


94 


12 27 20.69 


+23 19 47.5 


-14 


-10.0 


0.0 


ft 
u.u 


Q Ql. 

7.0J 


Q 46 


8 Q1 


8 S4 


8 4S 

O.'+J 


Q16 


64 


12 23 11.99 


+29 14 59.9 


-2.0 


-6.0 


0.0 


0.0 


10.29 


9.84 


9.13 


8.58 


8.47 


0.764 


0.79 


12 33 30.19 


+26 10 00.1 


\« 


-10.0 


0.0 



u.u 


in SQ 

iu. jy 


10 14 


Q 94 


8 79 

o. /z 


8 SQ 


766 
u. /oo 


7Q 

U. / 7 


12 39 52.43 


+25 46 33.0 


-zu.u 


-8.0 


0.0 


ft 
u.u 


1 O 7Q 
1U. ly 


1 O AQ 


Q 91 
y.zj 


8 74 
0. 


8 AS 


894 


64 


12 28 56.43 


+26 32 57.4 


-14 


-8.0 


0.0 





10 53 


10 25 


Q 91 

y.Zl 


8 77 
o. / / 


8 66 


852 


64 


12 24 53.60 


+23 43 04.9 


A ft 
-+.U 


-8.0 


0.0 


ft ft 
u.u 


i o 71 


1 O 14 

1U. J^t 


Q 10 

y. jy 


8 86 
O.OD 


8 89 
O.OZ 


O 840 
U.o'+U 


O 64 


12 28 34.29 


+23 32 30.6 


-8 


-14.0 


0.0 


ft ft 

u.u 


i o 40 

1U.H-U 


Q Q7 
y.y l 


9 47 


8 Ql 

O.y J 


8 86 


7Q8 
u. /yo 


79 


12 33 00.62 


+27 42 44.8 


-14 


-14.0 


0.0 





10 42 


9 80 


9 47 


8 94 


8 87 


805 


79 


12 35 17.03 


+26 03 21.8 


-2 


-10.0 


0.0 





10 69 


10 35 


Q S1 

y.J l 


Q 01 

7.U1 


8 Ql 

0.7J 


818 


64 


12 25 10.14 


+27 39 44.8 


-6.0 


-12.0 


0.0 


0.0 


10.69 


10.10 


9.57 


9.07 


8.93 


777 


0.79 


12 18 57.27 


+25 53 11.1 


-12 


-12.0 


3.0 


1 ft 

l.U 


1 O 80 
lU.oU 


1 O 10 


7.0H 


Q 08 


8 Q4 


798 
U. /Zo 


7Q 

U. / 7 


12 21 15.63 


+26 09 14.1 


-10.0 


-8.0 


0.0 


0.0 


10.87 


10.41 


9.62 


9.09 


8.97 


773 


0.79 


12 32 08.09 


+28 54 06.5 


-10.0 


-4.0 


0.0 


0.0 


10.74 


10.30 


9.62 


9.09 


8.99 


0.785 


0.79 


12 12 53.23 


+26 15 01.3 


-12 


-12.0 


0.0 





10 54 


Q 79 
y. I z 


9 58 


Q 1 1 

7.11 


8 99 


819 


64 


12 23 47.23 


+23 14 44.3 


-12 


-16.0 


0.0 





10 36 


9 54 


9 68 


9.13 


9 02 


759 


0.79 


12 18 17.77 


+23 38 32.8 


f, ft 
-u.u 


-14.0 


0.0 



u.u 


10 77 


10 28 


9 76 


9 20 


9 10 


7SQ 

U. i Jy 


7Q 

U. / 7 


12 22 52.37 


+26 38 24.2 


-8 


-10.0 


0.0 


0.0 


1 1 42 


1112 


9 78 


9 26 


9.11 


756 


79 


12 26 51.03 


+26 16 01.9 


-14 


-2.0 


0.0 





11 07 


1 n S6 


Q 


9 27 


9 16 


79 S 
U. /ZJ 


79 


12 09 12.44 


+26 39 38.9 


-16.0 


-6.0 


0.0 





10.60 


9.78 


9.83 


9.27 


9.18 


770 


0.79 


12 27 00.81 


+29 36 37.9 


-4.0 


-6.0 


0.0 


0.0 


11.07 


10.70 


9.80 


9.33 


9.20 


0.806 


0.79 


12 23 28.21 


+25 53 39.9 


in n 


-12.0 


1.0 


1 fi 

l.U 


1 1 43 


1 1 no 


Q 09 


Q IS 


Q 96 


7S8 
U. / Jo 


7Q 
U. / y 


12 24 10.37 


+29 29 19.6 


f\ 

-U.U 


-2.0 


0.0 


ft ft 
u.u 


i n 77 

1U. 1 1 


10 97 

1U.Z / 


10 06 
1U.UU 


Q SO 


Q 11 
y.J J 


681 

U.UO J 


7Q 

U. / 7 


12 34 46.93 


+24 09 37.7 


-12 


-4.0 


2.0 


5 ft 


1 1 18 

1 1 .JO 


1 O 87 

1U.0 / 


1 06 


Q SI 


Q 1Q 

y.Jy 


74Q 
U. I"ry 


7Q 

U. / 7 


12 15 34.01 


+26 15 42.9 


-12 


-6.0 


0.0 


ft ft 
u.u 


I 1 08 

I I .Wo 


1 O 96 


10 11 


S7 
y.J 1 


Q 47 

7.'4 / 


7S7 

U. / J / 


7Q 

U. / 7 


12 26 00.26 


+24 09 20.9 


-10 


-4.0 


2.0 


2.0 


13.85 


12.51 


10.98 


10.36 


10 14 


558 


80 


12 28 57.67 


+27 46 48.4 


-14 


-2.0 


3.0 


? ft 
z.u 


13 04 


1 1 37 


1 QQ 

1U.77 


1 IS 

1 U.J J 


10 19 


SS9 

U. J JZ 


80 


12 16 00.86 


+28 05 48.1 


-12 


-4.0 


2.0 


2 


14 15 


1 1 43 


11 07 


10 52 


10 24 


S41 


80 


12 30 57.39 


+22 46 15.2 


-12 


-4.0 


3.0 


1 ft 
1 .u 


14 17 


12 13 


1 1 24 


1 6S 


10 42 


SI 4 
U. J 1H- 


o so 

U.oU 


12 31 57.42 


+25 08 42.5 


- 1U.U 


-12.0 


0.0 


1 ft 
1 .u 


14 22 


12 94 


1 1 40 


1 7Q 
1U. ly 


10 SS 


4Q9 
U.H-yZ 


SO 

U.oU 


12 31 27.72 


+25 23 39.9 


R 
-o.u 


-6.0 


1.0 


? ft 
z.u 


1 4 07 


1 9 Q6 
IZ.yo 


1 1 44 


1 84 
iu.oh- 


1 61 

1U.UJ 


478 
U.H- / o 


SO 

U.OU 


12 25 55.76 


+29 07 38.3 


8 ft 
-O.U 


-18.0 


2.0 


1 ft 


11 81 


1 1 78 


1 1 S6 

1 1 .JO 


1 O QS 
lU.yj 


1 O 7S 

1U. / J 


O 460 
U.4-OU 


so 

U.oU 


12 23 55.53 


+23 24 52.3 


to n 


-4.0 


1.0 


ft ft 
u.u 


1 4 17 
i *+. J / 


12 64 


1 1 SQ 

1 1 .J7 


1 QQ 

1U.77 


1 77 

1U. / / 


4SS 


so 

U.OU 


12 25 02.64 


+26 42 38.4 


9. ft 
-O.U 


-4.0 


0.0 


^ ft 

J.U 


1 A 99 
14-. ZZ 


1 i ss 


1 1 69 
1 1 .OZ 


1 1 01 

1 l.U J 


1 O 7Q 
11/. ly 


O 4S9 


SO 
U.oU 


12 30 04.87 


+24 02 33.9 


- 1 


-8.0 


3.0 


ft ft 
u.u 


14 76 


13 27 


11 77 


1118 


10 94 


498 


80 


12 18 12.77 


+26 49 15.6 


-8 


0.0 


1.0 


1 


15.61 


12 78 


12 02 


1 1 46 


IMS 


0.392 


60 


12 15 16.93 


+28 44 50.0 


-12 


-14.0 


4.0 


1 


14 1 1 



u.u 


12 00 


1 1 35 


1 1 17 


18Q 
U.Joy 


60 

U.OU 


12 23 12.03 


+23 56 15.1 


-8.0 


-6.0 


2.0 


1.0 


15.40 


13 31 


12.20 


1161 


1 1 38 


0.352 


0.60 


12 31 00.28 


+26 56 25.1 


16 


-14.0 


3.0 


1 ft 
1 ,u 


14 04 


12 52 


12 25 


1 1 S6 
1 1 .JO 


1 1 1Q 

1 1 . jy 


1S1 
U.J J 1 


60 

U.OU 


12 33 31.35 


+24 12 09.1 


a n 


-16.0 


1.0 


? ft 
z.u 


14 76 

in-. /o 


1 1 4S 


12 29 


1 1 SQ 

1 1 J7 


1 1 40 


148 
U. jHo 


60 

U.OU 


12 16 37.30 


+26 53 58.2 


in n 


-6.0 


2.0 


1 ft 

J.U 


1 S 1 1 


12 94 


12 23 


1 1 68 

1 1 .UO 


1 1 42 


146 


60 

U.OU 


12 24 10.89 


+23 59 36.4 


f\ ft 

-O.U 


-4.0 


4.0 


1 ft 
1 ,u 


i s fin 

1 J.DU 


1 1 67 
i j.o / 


1 9 97 

iz.z / 


1 1 66 

1 l.OD 


1 1 4S 
1 1 .40 


O 140 
U. J'+U 


O 60 

U.OU 


12 28 38.70 


+25 59 13.0 


-6.0 


-2.0 


1.0 


2.0 


14.31 


13.94 


12.36 


11.69 


11.53 


0.326 


0.60 


12 16 22.84 


+24 19 01.1 


-12.0 


-2.0 


2.0 


4.0 


14.78 


13.96 


12.39 


11.73 


11.55 


0.323 


0.60 


12 27 08.56 


+27 01 22.9 


-16.0 


-2.0 


2.0 


3.0 


14.14 


12.52 


12.38 


11.73 


11.57 


0.319 


0.60 


12 28 04.54 


+24 21 07.6 


-12.0 


-10.0 


1.0 


3.0 


15.75 


14.29 


12.39 


11.84 


11.58 


0.318 


0.60 


12 38 04.72 


+25 51 18.5 


-16.0 


-6.0 


4.0 


0.0 


14.96 


13.43 


12.46 


11.84 


11.64 


0.308 


0.60 


12 14 19.78 


+25 10 46.6 


-6.0 


-12.0 


3.0 


3.0 


14.81 


13.92 


12.50 


11.80 


11.65 


0.305 


0.60 


12 28 50.08 


+27 17 41.7 


-20.0 


-12.0 


1.0 


2.0 


14.46 


12.84 


12.55 


11.84 


11.70 


0.297 


0.60 


12 36 34.30 


+25 00 38.3 


-14.0 


-2.0 


4.0 


6.0 


15.59 


13.43 


12.51 


11.96 


11.70 


0.298 


0.60 


12 26 37.32 


+22 34 53.4 


-12.0 


-8.0 


1.0 


1.0 


15.19 


13.34 


12.60 


11.97 


11.77 


0.288 


0.60 


12 33 30.31 


+28 12 55.9 


-12.0 


-14.0 


2.0 


11.0 


14.06 


13.59 


12.65 


12.07 


11.84 


0.279 


0.60 


12 16 29.21 


+23 32 32.9 


-12.0 


-8.0 


6.0 


1.0 


15.08 


14.55 


12.77 


12.13 


11.91 


0.270 


0.60 


12 30 46.17 


+23 45 49.0 


-10.0 


-6.0 


1.0 


9.0 


15.43 


13.83 


12.72 


12.17 


11.91 


0.269 


0.60 


12 19 37.99 


+26 34 44.7 


-6.0 


-8.0 


4.0 


2.0 


16.46 


14.03 


12.78 


12.24 


11.92 


0.269 


0.60 


12 14 23.97 


+28 21 16.6 


-4.0 


-8.0 


2.0 


1.0 


15.52 


0.0 


12.83 


12.14 


11.92 


0.269 


0.60 



6 S. L. Casewell et al. 



~i — 1 — i — 1 — i — 1 — i — 1 — i — 1 — r 
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1 0.2 0.3 0.4 
H-K 




K magnitude 



Figure 3. A JHK colour-colour diagram for candidate (crosses) and pre- 
viously known (solid triangles) cluster members. Empirical dwarf (dotted 
line) and giant (dashed line) sequences of Koornneef (1983), are over- 
plotted. The diagram has been plotted using the 2MASS system. The 
colour transformations of Carpenter (2001) were used to convert the model 
colours. 



panding the search radius to greater than 10 mas yr _1 unearthed no 
statistically significant increase in the number of candidate mem- 
bers, as the candidates and control stars increased proportionally, 
leading to many candidates with extremely small probabilities of 
being cluster members. Our survey is complete to 90% at this ra- 
dius, as explained earlier. As the stars get fainter however, the errors 
in their proper motion do increase - which is to be expected. It is 
entirely possible that the completeness is less than 90% for the faint 
stars. 

We have estimated the masses of our candidates and the previ- 
ously known cluster members using the evolutionary tracks of Gi- 
rardi et al. (2002), for masses^ 1 M , and the NEXTGEN models 
of Baraffe et al. (1998) for M<1 M . The stars were binned ac- 
cording to K magnitude or J-K colour, for the vertical and diagonal 
portions of the main sequence. We then linearly interpolated be- 
tween the model masses to estimate the masses of the cluster stars 
and our candidates. These masses are shown in the final columns of 
Tables 1 and 2, and illustrated in Figure 5. By multiplying the es- 
timated masses of the stars by their membership probabilities (we 
assume P m cmbcrship=l -0 for previously known members) and sum- 
ming, we determine a total cluster mass of ~102 M . This in turn 
allows us to derive a tidal radius of 6.5 pc or 4. 1° at 90 pc. Thus we 
find that within our adopted survey radius of 4° we should expect 
to find 99 per cent of the gravitationally bound proper motion se- 
lected cluster members. Indeed, increasing the search radius to 5°, 
led to a near equal increase in both candidate cluster members and 
stars in the control. 

Figure 3 shows the J-H, H-K two colour diagram for all new 
candidate and previously known cluster members. The colours of 
our candidates show a clear turnover at H-K Ri0.15, consistent with 



Figure 4. Luminosity function for the cluster, taking into account probabil- 
ity of membership. The error bars indicate the Poisson error. 




-0.1 0.1 

log mass 



0.2 0.3 



Figure 5. Mass function for the cluster, taking into account probability of 
membership. The error bars indicate the Poisson error. 



them being dwarfs and supporting our conclusion that these are low 
mass members of Melotte 111. Therefore, we have constructed a 
Ks luminosity function, and mass function for the cluster (Figures 
4, 5). The luminosty function is seen to peak at Ks =7.5, Mk s =2.7 
and to gradually decline towards fainter magnitudes. Superimposed 
on this decline is the Wielen dip (Widen 1971) at Ks=10.5. 
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The peaking of the luminosity function at this comparatively 
bright magnitude strongly suggests that a large proportion of the 
lower mass members, which are normally more numerous than the 
higher mass members (e.g. Hambly et al. 1991), have evaporated 
from the cluster. This should not come as a great surprise since at 
ages r<; 200Myrs, dynamical evolution of a typical open cluster 
is expected to lead to the loss of a significant and increasing frac- 
tion of the lowest mass members (e.g. de la Fuente Marcos & de 
la Fuente Marcos 2000). The mass function also supports this con- 
clusion. Thus despite having found low mass stars in the cluster 
our study does not contradict the conclusion of previous authors 
that the cluster is deficient in low mass stars. Odenkirchen et al. 
(1998) found that the luminosity function of extratidal stars associ- 
ated with the cluster rises towards lower masses. This can be con- 
sidered to be indicative of the ongoing loss of low mass members. 
Ford et al (2001) however, determined that ~ half of these stars 
did not belong to the cluster. A more in depth study of lower mass 
members would be required to determine a more accurate luminos- 
ity function at fainter magnitudes. 

Our faintest candidate member has an estimated mass of 0.269 
Mq. There are some 18 stars in the magnitude range Ks=ll-12, 
the lowest luminosity bin, thus it seems a distinct possibility that 
the Coma Berenices open cluster will have members right down to 
the limit of the hydrogen burning main sequence and possibly into 
the brown dwarf regime. This conclusion is also confimed by the 
rising mass function at lower masses. 



4 CONCLUSIONS 

We have performed a deep, wide area survey of the Coma Berenices 
open star cluster, using the USNO-B1.0 and the 2MASS Point 
Source catalogues, to search for new candidate low mass mem- 
bers. This has led to the identification of 60 objects with proba- 
bilities of cluster membership, P m embership^0.6. Our lowest mass 
new candidate member has M~0.269 Mq in contrast to the previ- 
ously known lowest mass member, M~0.915 M©. Thus we have 
extended considerably the cluster luminosity function towards the 
bottom of the main sequence. As reported by previous investiga- 
tions of Melotte 111, the luminosity function is observed to decline 
towards fainter magnitudes, indicating that the cluster has proba- 
bly lost and continues to loose its lowest mass members. This is 
not surprising for a cluster whose age is 400-500Myrs. Neverthe- 
less, as the cluster luminosity function remains well above zero at 
Ks=12, we believe the cluster probably has members down to the 
bottom of the hydrogen burning main sequence and possibly some 
brown dwarfs. Thus a deeper IZ survey of the cluster could prove a 
fruitful undertaking. 
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